Non-invasive vibrational mode spectroscopy of ion Coulomb crystals through resonant 

collective coupling to an optical cavity field 
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We report on a novel non-invasive method to determine the normal mode frequencies of ion 
Coulomb crystals in traps based on the resonance enhanced collective coupling between the electronic 
states of the ions and an optical cavity field at the single photon level. Excitations of the normal 
modes are observed through a Doppler broadening of the resonance. An excellent agreement with the 
predictions of a zero-temperature uniformly charged liquid plasma model is found. The technique 
opens up for investigations of the heating and damping of cold plasma modes, as well as the coupling 
between them. 
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Cold one-component plasma physics [l[ has in the past 
two decades led to a series of interesting results due to 
the availability of fast computers as well as the pos- 
sibility to experiment with ensembles of trapped, laser- 
cooled atomic ions Prominent examples are the 



understanding of structural properties of crystallized cold 
plasmas in both Penning @, Q and Paul [Mil [H EH 
traps, and the investigation of the normal mode dynam- 
ics of cold magnetized plasmas in Penning traps llHl4|. 



While there exist many similarities between experiments 
in Penning and Paul traps, the unmagnetized plasmas in 
Paul traps are e.g. known to heat up much faster than the 
magnetized plasmas in Penning traps due to the presence 
of the rf fields [1 , 17 1 . The lack of a rotational symmetry 
axis in Paul traps has as well been found to be responsible 
for the observation of specific crystalline structures fljj . 
Exploration of the normal mode dynamics of unmagne- 
tized ion Coulomb crystals in e.g. linear Paul traps will 
hence add to our understanding of the influence of the 
trapping environment on the physics of such crystals. 
Furthermore, since large Coulomb crystals are excellent 
candidates for the realization of high-fidelity quantum 
memories for light [l8], such studies can shed light on 
the influence of the excitation of these modes on the fi- 
delity as well as on the prospect of storing several pho- 
tonic quantum bits through deliberate excitations of spe- 
cific vibrational modes. Larger Coulomb crystals have as 
well recently been considered as a system for performing 
quantum simulations, in which respect knowledge of nor- 
mal mode dynamics is needed (l9| . Finally, ion Coulomb 
crystals represent extremely interesting systems to study 
cavity optomechanics phenomena [20l. l2l|. since, in spite 
of their solid nature, they possess free atomic resonance 
properties and can hence be made very sensitive to the 
radiation pressure force exerted by optical fields. 

In this Letter, we report on the study of normal mode 
vibrations of Coulomb crystals of 40 Ca + ions in a lin- 
ear Paul trap by a novel non-invasive technique. The 
technique is based on monitoring the coherent collec- 



tive resonant response of the atomic ions constituting 
the crystal to a single photon optical cavity field fljj . 
By having a standing wave optical cavity incorporated 
in the trap setup with the light propagation axis coincid- 
ing with the rf-field free axis of the trap, we can detect 
very small changes in the optical response due to Doppler 
shifts resulting from the ions' motion along the cavity 
axis when a mode is excited. Since the probing does not 
rely on the observation of incoherently scattered pho- 
tons as e.g. the case in Doppler velocimetry ll-lj| and 
in Sympathetically-Cooled Single Ion Mass Spectrome- 
try 22j, the detected signal can in principle be purely 
dispersive in nature and does not require any excitation 
of the ions. The technique is not limited to linear Paul 
traps, but should work as well in any open structured 
trap geometries such as e.g. the Penning-Malmberg type 
traps [Tl|, [23]. 

Previous Doppler velocimetry imaging experiments 
with cold magnetized plasmas of Be + ions in a Penning 
trap have led to the identification of a series of normal 
modes [ll]4l3 |. which correlate to specific (7,m)-modes 
theoretically predicted for spheroidal shaped uniformly 
charged liquids 24]. For the cold unmagnetized plasmas 
in the linear Paul trap used in our investigations, the 
corresponding low-order normal mode frequencies are ex- 
pected to be close to those of the (I, m) modes, which are 
given by 24] 
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where lu p = yj q 2 no/Meo is the plasma frequency (q is the 
charge of the ion, M its mass and no the uniform ion den- 
sity), If 1 = P ; m (l/Vl -a 2 ) and Qf = Qf(l/Vl -a 2 ) 
are first and second-order Legendre polynomials with 
cylindrical indices (I, m) and with the prime denoting dif- 
ferentiation with respect to the entire argument. Finally, 
a is the aspect ratio of the plasma given by a = R/L, 
where L and R are the half-length and radius of the 
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FIG. 1: (color online.) Schematic of the linear Paul trap with 
an integrated optical cavity. The excitation of the crystal 
modes (AC modulation) can be monitored either by measur- 
ing the change in the reflection spectrum of a weak cavity field 
(866 nm) with an avalanche photodiode (APD) or by imag- 
ing the fluorescent light (397 nm) emitted by the 40 Ca + ions 
during laser cooling onto a CCD camera. The inset shows the 
relevant energy levels and transitions of the 40 Ca + ion. 



spheroid, respectively (see also Fig. [2^,). 

The corresponding spatial structure of the modes can 
generally be expressed in terms of generalized spheroidal 
coordinates [25[ . While the modes generally have a non- 
trivial dependence on the ions' position in the crystal, 
for longitudinal modes (l,m = 0), the displacement 5z, 
close (r ~ 0) to the axis of rotational symmetry (z-axis), 
is approximately given by 



Sz oc P?'(z/y/L 2 -R 2 ). 



(2) 



dant calcium 2(| 28 1. The Coulomb crystals are created 
through Doppler laser-cooling along the trap axis by two 
counter-propagating beams at 397 nm (total power of 
~ 4 mW and beam diameter of ~ 1 mm) and an 866 
nm beam applied from the side to prevent the ions from 
being shelved into the D 3 / 2 state (see insert of Fig. [TJ). 
The number of ions constituting the crystals can be de- 
duced by imaging the light (397 nm) scattered by the ions 
during the cooling process onto a charge coupled device 
(CCD) chip [H. 

Normal modes of the ion Coulomb crystals can be ex- 
cited by applying identical AC electrical potentials at 
variable frequency to the four end-pieces at each end of 
the trap. By this simple geometry modes corresponding 
to the (I, m = 0) charged liquid modes can easily be ex- 
cited by having the voltages applied at each end either 
in-phase (I even) or 180° out of phase (I odd). 

To monitor the coherent collective resonant response 
of the atomic ions constituting the crystal, a 11.8 nun- 
long optical cavity with a measured TEMoo mode waist 
of 37 /im is used [2f|. The optical finesse of the cavity 
is about 3,000 at the 866nm resonant wavelength of the 
3d 2 D 3/2 - 4p 2 P 1/2 transition in the 40 Ca+ ion addressed 
in the experiments. 

The plasma mode diagnostic involves a sequence of 
steps. First, a Doppler cooling period of 5 /is is fol- 
lowed by a period of 12 /is, where the ions in the crystal 
are prepared in the mj = +3/2 magnetic sub-state of 
the long-lived metastable D 3 / 2 level by optical pumping 
(efficiency ~ 97%). Next, a weak left-handed circularly- 
polarized pulse (1.4 /is long) of 866 nm light is coupled 
into the cavity to probe the collective response of the ions. 
The mean intracavity photon number is less or about one 
at any time. During this probing period, the photons re- 
flected by the cavity are measured by an avalanche pho- 
todiode with an overall collection efficiency of 16%. This 
sequence is repeated at a rate of 50 kHz while the cavity 
length is scanned at a rate of 30 Hz. The cavity lineshape 
is reconstructed by averaging a few hundred scans (see 
fl8j for more details). 

With the ion Coulomb crystal at rest, the cavity probe 
linewidth is given by fl8j ]: 



Consequently, all modes with m = have a spatial vari- 
ation along the z-axis which, under excitation, will lead 
to measurable Dopplcr-shifts of the ions resonance fre- 
quency along this axis. 

Figure [1] shows a schematic of the experimental setup 
for our studies (a more detailed description of the cavity- 
trap setup can be found in Ref. t 26]). The linear Paul 
trap consists of four cylindrical electrodes to which an 
rf voltage U r f between 200 V and 300 V is applied at 
a frequency ft of 27rx4 MHz for radial confinement (xy- 
plane in Fig. [TJ. To obtain a linear quadrupole config- 
uration, the rf voltages are 180° out of phase between 
nearest neighbor electrodes. Each of the four electrodes 
are sectioned into three parts to allow for the application 
of DC (1-10 V) and AC voltages to the eight end-pieces 
for static confinement of the ions along the cavity axis 
(z-axis) and for normal mode excitations, respectively. 
The electrodes have diameters of 5.2 mm and the near- 
est neighbour electrode spacing is 1.8 mm. With these 
parameters, the single ion effective (pseudo) potential is 
to a very good approximation harmonic in all dimensions 
in the region of the crystals, which leads to constant den- 
sities of the crystallized 40 Ca + ions with values between 
2.6 x 10 8 and 5.7 x 10 8 cm" 3 [27|. 

The 40 Ca + ions are produced through isotope selec- 
tive photoionization of atoms in a beam of natural abun- where k is the wavevector of the cavity field. 
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where k is the cavity decay rate, gjq the collective cou- 
pling rate, 7 the optical dipole decay rate and A the 
detuning of the probe with respect to atomic resonance. 
If the crystal is moving with a fixed velocity v z along the 
cavity axis the probe linewidth would due to the Doppler 
shift be given by 



k' = k + g 2 N 



7 (7 2 + A 2 + (£^) 2 ) 
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(4) 



3 




139 140 141 142 143 144 145 

excitation frequency (2n kHz) 



FIG. 2: (a) Projection image of a 1.2 mm-long crystal with 
4000 ions and density 2.6 x 10 8 cm" 3 . The DC and rf voltage 
amplitudes are 2.36 V and 205 V, respectively, and the aspect 
ratio q = R/L is 0.135. (b) Cavity probe linewidth as a 
function of the mode excitation frequency applied to drive 
the equivalent of the (2, 0) mode of this crystal. The probe is 
tuned to atomic resonance (A = 0) and the AC modulation 
voltage is 0.75 V. The red line is a Lorentzian fit. 

When exciting one of the normal modes, all the ions 
in the crystal will generally not have the same velocity 
at a given instance (this is only true for the center-of- 
mass mode, see Eq. ((2)), and, furthermore, the velocity 
of the individual ions will change during the mode-period. 
However, if the timescale at which the velocity of the in- 
dividual ions changes is slow as compared to the effective 
cavity photon decay time, one can model the expected 
probe lineshape by averaging the contributions from the 
different parts of the crystal over a single mode-oscillation 
period. From Eqs. ((3|) and (j4j it follows that any excita- 
tion of the ions' motion will generally lead to a narrower 
linewidth due to the reduced coupling strength. 

Fig. (2b shows the cavity probe linewidth for A=0 
as a function of the modulation frequency applied to 
excite the equivalent of the (2,0) "quadrupole" mode 
of a 1.2 mm-long crystal with 4000 ions at a density 
2.6 x 10 8 cm -3 (Fig. [2^). A clearly reduced linewidth is 
observed around 142 kHz. The exact resonance value of 
w^ e o y/27r = 142.1 ± 0.1 kHz is in good agreement with 

the resonance frequency Lu^°^ el /2tt — 142.2 ±1.1 kHz 
expected for the (2, 0) mode of a charged liquid crys- 
tal for this aspect ratio and charge density. In order to 
test more generally how well the uniformly charged liquid 
model describes unmagnetized ion plasmas confined in a 
linear Paul trap, the resonance frequencies of the lowest 
lying normal modes of ion crystals with various aspect 




aspect ratio a 



FIG. 3: Resonance frequencies corresponding to the (1,0) 
(squares), the (2,0) (circles) and the (3,0) (triangles) modes 
as a function of the aspect ratio a, for a fixed plasma fre- 
quency uj p /2n — 536 kHz (U r f = 205V^). The solid lines show 
the theoretical predictions of Eq. (fljl. The error bars are 
within the point size. The inset shows the (2, 0) mode reso- 
nance frequency as a function of the (2, 2) mode off-resonant 
modulation depth £[/ 2 ) ( see text) for a fixed aspect ratio 
a — 0.135. The mode resonance frequency is normalized 
to that expected for a plasma without any excitation of the 
(2, 2) mode and the red lines show the uncertainty in the 
expected resonance frequency due to the density calibration. 
The Coulomb crystals contain between 4000 and 12000 ions, 
with densities of 2.6 — 5.6 x 10 8 cm -3 and their temperature 
is in the 10 mK range. 



ratios have been determined by monitoring the linewidth 
of the cavity probe field when tuned to the atomic res- 
onance. The resulting mode resonance frequencies are 
presented in Fig. [3] together with the predicted values 
from the charged liquid model (see Eq. ([T])). The mea- 
sured values are consistent with the theory to better than 
one percent for all experimental data. This accuracy may 
appear quite surprising considering that during all these 
(1, 0) mode measurements, the (2, 2) mode is continu- 
ously off-resonantly excited by the linear rf quadrupole 
field confining the plasma. The rf-induced modulation 
depth is indeed up up to 20% of the radial extension 
of the crystal, i.e. comparable to that typically used 
for the axial excitation of the (Z,0) modes. The inset 
of Fig. [3] clearly shows, however, that, within the cur- 
rent experimental accuracy, the (2,0) mode frequency 
does not have any systematic dependence on the (2, 2) 
mode modulation depth £^/ 2 )' defined by the micromo- 
tion amplitude of an ion at the radial position (xo,yo) 
(x(t),y(t)) = (x (l + C { 1 2 ) cos fit), y (l - cos^)). 
This result is though consistent with molecular dynam- 
ics simulations from which it has been predicted that the 
radial rf field-driven micromotion in linear Paul traps 
should have an extremely weak coupling into the axial 
motion of the ions [3f . 
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FIG. 4: Cavity probe linewidth as a function of the probe 
detuning, with (circles) or without (squares) modulation at 
the (1,0) mode resonance frequency, for a crystal similar to 
the one presented in Fig. O The solid lines are fits according 
to the theoretical model described in the text. 



In Fig. 21 the cavity probe linewidth with and with- 
out exciting the (1,0) mode at the resonance frequency 
o)/27r = 94 kHz is presented as a function of the 
detuning of the probe with respect to the atomic res- 
onance. From the Lorentzian line profile in absence of 
the mode excitation, one finds (k, 7, gjv)/27r = (2.2 ± 
0.1, 12.3 ± 0.3, 8.2 ± 0.1) MHz. When the (1, 0) mode is 
excited, the absorption line is substantially modified to 
a non-Lorenzian profile. A fit based on the (1,0) mode 
function (Eq. (j4])) averaged over one oscillation period 
results in a driven motion amplitude v jf = 5.9 ± 0.7 
m/s (v z (t) = v z o cos(w( 10 )i + f))- This value is in good 



agreement with the value v 



fluo 
zO 



5.3 ± 0.6 m/s de- 



duced from fluorescence images recorded in phase with 
the modulation voltage. Here, a motional amplitude 
z fiuo _ g _|_ ^ ^ m wag measure( j a t the mode resonance 

frequency LJnm/2'TT — 94 kHz. Since equally good quanti- 
tative agreement between experiments and the model was 
also found for the (2, 0) mode, this proves that quantita- 
tive information about the ions' motion can reliably be 
obtained from the ion-cavity coupling without the need 
for observing directly the fluorescence signal. In the fu- 
ture, higher order (/, 0) modes are envisioned to be stud- 
ied through mode excitations using spatial- and time- 
modulated radiation pressure forces. For 40 Ca + ions this 
can e.g. be achieved through the combined application 
of a 866 nm repumper beam with a spatially modulated 
intensity profile along the cavity axis and a time- varying 
intensity of one of the 397 nm cooling beams 22[ . Even- 
tually, for the high spatial modulation of modes with 
large I, the liquid model should cease to apply. Further 
applications could be measurements of ion Coulomb crys- 
tal temperatures and heating rates [14| and more detailed 
investigations of the coupling between the various normal 
modes at various temperatures and structural phases of 



the plasma d[25|. Finally, we also believe that the spec- 
troscopic findings as well as the non-invasive character of 
the method used will be important for e.g. obser ving ra- 
diation pressure-induced optomechanical effects 20J, |l 
with ion Coulomb crystals. 
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